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ABSTRACT 

The TASTE project is searching for low-mass planets with the Transit Timing Variation (TTV) technique, by gathering high-precision, 
short-cadence light curves for a selected sample of transiting exoplanets. It has been claimed that the "hot lupiter" WASP-3b could 
be perturbed by a second planet. Presenting eleven new light curves (secured at the IAC80 and UDEM telescopes) and re-analyzing 
thirty-eight archival light curves in a homogeneous way, we show that new data do not confirm the previously claimed TTV signal. 
However, we bring evidence that measurements are not consistent with a constant orbital period, though no significant periodicity can 
be detected. Additional dynamical modeling and follow-up observations are planned to constrain the properties of the perturber or to 
put upper limits to it. We provide a refined ephemeris for WASP-3b and improved orbital/physical parameters. A contact eclipsing 
binary, serendipitously discovered among field stars, is reported here for the first time. 

Key words, techniques: photometric - stars: planetary systems ~ stars: individual: WASP-3, VSX J183407. 3+353859 



1. Introduction 

Most of the extrasolar planets discovered so far are massive, 
gaseous giant planets. The present trend is to probe smaller and 
smaller masses, with the final aim of detecting temperate "super- 
Earths" or Earth-sized rocky planets (1-10 Me, 1-3 7?®) around 
solar- type stars. The signal expected from a true Earth analog 
orbiting a Sun twin is extremely small: ~ 80 ppm for the photo- 
metric transit, and ~ 10 cm/s for the radial velocity (RV) Doppler 
shift. In the conventional framework both measurements are re- 
quired to derive the planetary radius and mass (Rp, Mp), i.e. the 
basic quantities necessary to confirm the planetary status of the 
transiting body, and to characterize it. Achieving a photometric 
precision of ~ 20 ppm (which would allow the detection of an 
Earth-like planet at 4cr level) over the timescale of a transit is 
within reach of space-based telescopes only, while a long-term 
< 20 cm/s RV accuracy is still too ambitious even for the most 
stable spectrographs, like HARPS. 

A few indirect techniques have been developed to get es- 
timate of Mp (or upper limits to it) for a non-transiting planet 
without the need of RV measurements, the most promising be- 
ing the Transit Timing Variation (TTV) method. By monitoring a 
known transiting planet by high-precision photometry, the cen- 
tral instant To of each individual transit can be estimated. The 
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gravitational perturbation of a previously-unknown third body, 
not necessarily transiting, can cause a significant variation of 
the orbital period P (Holman & Murray 2005). The effect is 
greatly increased if the perturber is locked in a low-order mean- 
motion resonance with the transiting planet (Agol et al. 2005). 
The Asiago Search for Transit timing variation of Exoplanets 
(TASTE) project was started in 2010 to search for TTVs with 
several ground-based, medium-class facilities on a sample of 
carefully selected targets (Nascimbeni et al. 201 la). 

The TTV technique has been exploited and already gave 
interesting results on a number of Kepler mission candidates. 
Some candidate planets in multiple systems were validated 
through TTV analysis (Lissauer et al. 2011). In the past few 
years, some authors have claimed TTV detections also from 
ground-based facilities, but none is confirmed so far. Among 
these claims, WASP-lOb (Maciejewski et al. 2011), WASP-5b 
(Fukui et al. 201 1), the intriguing case of HAT-P-13b, that was 
monitored also by TASTE (Pal et al. 2011; Nascimbeni et al. 
2011b; Southworth et al. 2012), and the subject of the present 
investigation: WASP-3b (Maciejewski et al. 2010). 

WASP-3b is a typical short period (P ^ 1.8468 days) 
"hot Jupiter" (1.31 /?jup, 1.76 Mjup), hosted by a F7-8 dwarf. It 
was discovered by Pollacco et al. (2008). Analyzing an (9 - C 
{Observed - Calculated) diagram computed by comparing the 
T() of fourteen transits with the value predicted by a linear 
ephemeris, Maciejewski et al. (2010) claimed the detection of 
a sinusoidal modulation with a period of Pjtv - 127.4 days and 
a semi-amplitude of ~ 0.0014 days ^ 2 min. They interpreted 
this signal as the effect of an outer perturber, identifying three 
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possible orbital solutions in the range 6-15 and P = 3.03- 
3.78 days. No independent confirmation of this claim has been 
published so far, though both Littlefield (2011) and Sada et al. 
(2012) discussed the consistence of their data with that TTV 
modulation. 

In this paper we present (Section 2) eleven unpublished tran- 
sits of WASP-3b: six of them have been gathered at the lAC- 
80 telescope and five at the UDEM 0.36m (Universidad de 
Monterrey, Mexico). We also sifted the literature in search of all 
the archival photometric data useful for a TTV study (Section 
3). In Sections 4 and 5 we describe how both new and archival 
light curves, for a total of forty-nine transits, were reduced and 
analyzed in a homogeneous way, with the same software tools 
to provide a consistent estimate of the planetary parameters and 
their uncertainties. This is crucial especially for To, whose es- 
timate has been shown to be easily biased by the employed 
analysis technique (Fulton et al. 2011; Southworth et al. 2012). 
Besides To, we also refined the orbital and physical parameters 
of WASP-3b, and computed an updated ephemeris (Eq. 10) for 
any forthcoming study on this target. In Section 5.2 we demon- 
strate that the TTV claimed by Maciejewski et al. (2010) is not 
supported by our analysis, and probably due to small-sample 
statistics. Yet, we point out that the revised O - C diagram dis- 
plays a complex, non-periodic structure and is not compatible 
with a constant orbital period. Finally, in Section 6 we discuss 
the possible origin of this TTV signal, and show that careful dy- 
namical modeling and additional photometric and RV follow-up 
is required to confirm the hypothesis and to constrain the mass 
and period of the possible perturber(s). 

2. TASTE observations 

2.1. IAC-80 observations 

We observed six transits of WASP-3b between 2011 May 7 
and Aug 2, employing the CAMELOT camera mounted on the 
1AC80 telescope. A log summarizing dates and other quantities 
of interest is reported in Table 1. Individual transits are iden- 
tified by an ID code ranging from Nl to N6. IAC80 is a 0.8m 
Cassegrain reflector installed at the Teide Observatory (Tenerife, 
Canary Islands) and operated by Instituto de Astrofisica de 
Canarias (lAC). CAMELOT is a conventional imager with a 
10' .4 X 10' .4 field of view (FOV), equipped with a E2V 42- 
40 2048x2048 CCD detector, corresponding to a 0.304" pixel 
scale. The software clock interrogated to save the timestamps in 
the image headers is automatically synchronized with the GPS 
time signal. 

All the observations were carried out with a standard Bessel 
R filter and the same instrumental setup. Windowing and 2x2 
binning were employed to increase the duty-cycle of the pho- 
tometric series, as described in Nascimbeni et al. (2011a). A 
10.4' X 3.2' read-out window was chosen to image the target 
and a previously selected set of reference stars in a region of the 
detector free from cosmetic defects. The one-amplifier readout 
was preferred to prevent gain offsets between different channels. 
Stars were intentionally defocused to a FWHM of 10-13 binned 
pixels (^ 6.0-8.0") to avoid saturation and to minimize system- 
atic errors due to intra-pixel and pixel-to-pixel inhomogeneities 
(Southworth et al. 2009). Exposure time was set to 20 s (Nl-4) or 
15 s (N5-6), resulting in a net cadence t = 19-25 s and a ~ 75- 
87% duty-cycle, with the only exception of N3. On that night, 
due to software problems, the images were read unbirmed and 
in unwindowed readout mode, decreasing both signal-to-noise 
ratio (S/N) and duty-cycle. 



The weather was photometric on all nights, except for a few 
thin veils during N3. Our initial goal was to start the series one 
hour before the first contact of the transit, and to stop one hour 
after the last contact. Nevertheless, the amount of pre-transit 
photometry is only a few minutes for the N4 and N6 light curves, 
and the N 1 series was interrupted twenty minutes earUer because 
of twihght. 

2.2. UDEM observations 

We observed five transits of WASP-3b between 2008 Sep 3 
and 2010 Aug 15 with the Universidad de Monterrey (UDEM) 
0.36-m reflector. UDEM is a small private college observa- 
tory having Minor Planet Center Code 720 located in the sub- 
urbs of Monterrey, Mexico. The data were acquired using stan- 
dard Bessel /- and Sloan z-band filters with a 1280x1024 pixel 
CCD camera at 1 .0" pixel scale, resulting in a field-of-view of 
~ 21'. 3 X 17'. 1. The observations were slightly defocused to im- 
prove the photometric precision and to avoid saturation. On-axis 
guiding was used to maintain pointing stability. Exposure times 
were set to 30 s for / and 40 s for z- All images were binned 2x2 
to facilitate rapid readout (~ 3 s). Each observing session lasted 
about 4.5 hours in order to accommodate the transit event and 
also to cover about one hour before ingress and one hour after 
egress. The computer clock was reset to UTC via Internet at the 
beginning of every observing session to the nearest second. 



3. Archival light curves 

All the available archival data were searched for any light curve 
useful for a TTV analysis, i.e. complete transits of WASP-3b 
with a suitable S/N. Partial transits were rejected, except for a 
few cases with high S/N and without obvious presence of "red 
noise" after a visual inspection. We define red noise, foUowing 
Pont et al. (2006), as correlated noise having covariance between 
data points on time-scales of the same order of the duration of 
our transit signal. The transit itself must be modeled and re- 
moved from a Ught curve in order to study the statistical proper- 
ties of its noise content. Thus, first we analyzed the whole sam- 
ple of light curves in Section 5.1, then we discuss how to select 
a high-quality subset for furter analysis in Section 5.2 ("ALL" 
vs. "SELECTED"). 

Hereafter we refer to each Ught curve with the corresponding 
alphanumeric ID code reported in Table 1 . Overall, thirty-eight 
archival light curves were (re-)analyzed in the present study: 

- Gibson et al. (2008): two transits in 2008 (Gl-2) observed 
with the RISE camera mounted at the 2.0 m Liverpool 
Telescope. Both nights were photometric. The second series 
(G2) ended just before the last contact. The RISE camera fil- 
ter has a custom wideband 500-700nm, approximately cov- 
ering both the Johnson R and V passbands. 

- Tripathi et al. (2010): five transits in 2008/2009 (Tl-5) 
from two diflJ'erent telescopes (Fred Lawrence Whipple 
Observatory FLWO-1.2m and University of Hawaii UH- 
2.2m) and in three Sloan passbands (g, i, z), gathered to com- 
plement their high-precision RV measurements. The first two 
FLWO ; transits (Tl-2) are partial, the first one missing the 
first contact by only few minutes. 

- Damasso et al. (2010): one complete transit observed in 2009 
(Dl), with a coimnercial 0.25 m reflector in a nonstandard/? 
band, as a part of the feasibility study for a project dedicated 
to the search for transiting rocky planets around M dwarfs. 
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Fig. 1. Light curves of WASP-3b analyzed in this Paper. The ID code of each transit matches the corresponding entry in Table 1. Data points are 
plotted with the original cadence, except for Gl, G2 and Zl, E7 that are binned respectively on 30 s and 300 s intervals for clarity. The red line is 
the best-fit model found by JKTEBOP. Transits have been ofi"set in magnitude by integer multiples of 0.025. 
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- Maciejewski et al. (2010): six complete transits in the R band 
from two different telescopes (Ml-3 at Rohzen-0.6m, M4-6 
at Jena-0.9m). The first five observations (Ml-5) were car- 
ried out on focused images, a practice resulting in a possible 
increase of the content of red noise. 

- Christiansen et al. (2011): six complete (Cl-4, C7-8) and 
two partial (C5-6) transits, extracted from a 18-day photo- 
metric series secured by the High Resolution Imager (HRI) 
mounted on the spacecraft EPOXI, as part of the EPOCh 
project. No filter was employed. 

- Littlefield (2011): three transits observed with a 11" 
Schmidt-Cassegrain reflector and without filter (LI -3). 
These three light curves correspond to the first, third, and 
fifth one presented in the original paper. The S/N of the other 
curves is too low for the inclusion in the present study. 

- Zhang et al. (201 1): one 2009 transit in the V band from the 
Weihai-l.Om telescope (Zl). 

- Sada et al. (2012): four transits secured in 2009-201 1 by dif- 
ferent telescopes: SI and S3-4 at the 0.5m Visitor Center 
telescope at Kitt Peak National Observatory (KPNO), in the 
Sloan z band; S2 at the KPNO 2.1m reflector in the J band, 
exploiting flie FLAMINGO infrared imager. 

Other six light curves from various amateur observers 
(El -7) were visually selected and downloaded from the ETD 
database'. 

The first two follow-up light curves of WASP-3b were pub- 
lished by Pollacco et al. (2008): one transit observed at the 
IAC80 (alternating V and / filters) and one at the Keele-0.6m 
{R band), both in 2007. Unfortunately, the original light curves 
are no more retrievable and we did not include them in this study. 
A To data point is determined by the first term of their published 
ephemeris: 

ro(BJDuTc) = 2454143.8503!!!:!™^ + N ■ 1.846834!!!;!!S (1) 

One should be careful, however, as this Tq (corresponding to 
2007 Feb 12) is not an independent measurement of a single 
transit event. Instead, it comes from an ensemble analysis of 
SuperWASP, IAC80 and Keele-0.6m data. For this reason we 
adopted Tq from Eq. (1) as our first data point in the determi- 
nation of our new ephemeris, but not in our subsequent TTV 
analysis (Section 5.2). The epoch A'^ of each reported observa- 
tion in our set is defined according to the ephemeris in Eq. (1), 
i.e. as the number of transits elapsed since 2007 Feb 12. 

All timestamps were converted to BJD(TDB), i e. based 
on Barycentric Dynamical Time following the prescription by 
Eastman et al. (2010). Each light curve was calibrated in time 
by identifying the time standard reported in the data headers 
and (when clarifications were necessary) contacting the authors. 
Generally speaking, even when the authors report a reliable syn- 
chronization source (GPS, NTP, etc.) for their data, it is impos- 
sible to carry out an external confirmation of that. The only ex- 
ception, discussed in Section 6, is when two or more independent 
observations of the same event were performed. Though we did 
not find any reason to doubt about the accuracy of the absolute 
time calibration of the above-mentioned data, it is worth noting 
that most of these observations were not performed with the spe- 
cific goal of a TTV analysis. Thus, the precision achievable on 
To could be limited by sub-optimal choices about the instrumen- 
tal setup. For instance, it is widely known that a "large" (in our 
case, T > 60 s) exposure time is one of those limiting factors 
(Kipping 2010). 

' http : //var2 . astro . cz/ETD (Poddany et al. 2010) 



When required, light curves in flux units were converted to 
magnitudes and normalized to zero by fitting a low-order poly- 
nomial function to the off-transit data points. In a few cases, 
the tabulated photometric errors are underestimated up to about 
50%, as confirmed by the reduced Xr^^ix^ being defined as 
thex^ divided by the number of degrees of freedom of the fit). 
This is not unusual in high-precision photometry, due to effects 
which are not accounted for by standard noise models (Howell 
2006): poorly-modeled scintillation, high-frequency systematics 
mimicking random errors, stellar microvariability for both tar- 
get and comparison stars, etc. We dealt with this by rescaling 
the errors by a factor of s/x^, following a common practice 
(Winn et al. 2007; Gibson et al. 2009). When flie photometric 
errors were not pubhshed, the error was assumed to be constant 
and equal to the scatter cr of the ofl'-transit polynomial-corrected 
curve. We define the scatter as the 68.27* percentile of the resid- 
ual distribution from the median, after a 5(t iterative clipping. 
This measurement is much more robust against outliers than the 
classical RMS. 



4. Data reduction 

4.1. IAC80 photometry 

The six 1AC80 light curves (Nl-6) were reduced with the 
STARSKY photometric pipehne, presented in Nascimbeni et al. 
(2011a,b) but here upgraded with some improvements. The 
present version v 1.1.002 adopts a new, fully empirical weight- 
ing scheme to carry out differential photometry. Reference stars 
were previously weighted by the amount of scatter measured 
on their light curves after being registered to the total reference 
magnitude m, (that is, the weighted mean of the instrumental 
magnitudes of all the comparison stars; Broeg et al. 2005). That 
was an iterative process. Instead, the updated version first ex- 
tracts the off-transit part of the series, then constructs a set of 
"intermediate" light curves of the reference stars by subtracting 
the magnitude of each of them to the off-transit magnitude of 
the target. IdeaUy, those curves should be flat and their RMS 

should be equal to the quadratic sum .^d-^ + cr^, being cr, the 

theoretical photometric noise expected on the target (calculated 
as in Nascimbeni et al. 2011a) and cr the intrinsic noise of the 
reference star, defined as at the end of Section 3. We therefore 

estimated the latter as o" = ^cr^ - cr^. The individual weights 

for a given comparison star are then assumed to be 1 / Vd^. The 
output of this weighting algorithm is checked during each run 
against two other weighting schemes: 1) equal weights, i.e. un- 
weighted, and 2) using weights derived from the expected the- 
oretical noise computed for the reference stars. The DSYS and 
PSYS parameters (as defined in Nascimbeni et al. 2011b) allow 
us to diagnose "bad" reference stars, and to set their weights to 
zero. For all the N1-N6 transits, the same set of eleven reference 
stars was employed, for consistency reasons. All of them show 
no sign of variability or higher-than-expected scatter. 

A second improvement to STARSKY is a new algorithm de- 
veloped to deal with Ught curves having red noise caused by 
veils, trails or thin clouds. This happen when the cloud possess a 
structure at angular scales on the same order of the FOV, and/or 
it is moving fast. Even differential photometry can be affected by 
these events, as the change of transparency can afl'ect the target 
and the reference stars by a different amount. Of course, this sys- 
tematic effect is correlated on with the rapidity of transparency 
changes, and this correlation can be exploited to discard the af- 
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fected frames. A quantity that we called "numerical derivative 
of the absolute flux" (NDAF) is evaluated for each frame i of the 
series: 

NDAF=l ^'""'^"-^> (2) 

being F, the weighted instrumental flux of the reference stars 
(Fi = 10""''"'', m, defined as above), and f, the JD time at the 
mid-exposure. When NDAF deviates more than 4cr from its av- 
erage along the series, the frame is discarded. The first and last 
frame are ignored by our algorithm. 

STARSKY outputs light curves extracted from a set of dif- 
ferent photometric apertures. These curves are then detrended by 
a routine which searches for Hnear and polynomial correlations 
between the ofl'-transit flux and several combinations of external 
parameters such as: the position of the star on the detector, the 
airmass, the FWHM of the stellar profiles, the mean sky level, 
the reference flux F„ and time. In the following, we describe the 
reason to search for a linear correlation between differential flux 
and airmass X. Let's consider the simple case of two stars (tar- 
get "t" and reference "r") having out-of-almosphere magnitudes 
mi and and color and ^r, respectively. Given the extinction 
coefficient k and the extinction color term k', the measured dif- 
ferential magnitude Am is the difference between the observed 
magnitudes: 

Am = (mt + kX + k'X^t) - + kX + k'X^^) (3) 

By grouping the involved terms, and taking into account that X 
is the only quantity that shows short-term variations: 

Am = (m, - m,) + X-k'(£r- f,-) (4) 

It is easy to see that systematic efi'ects due to diff'erential extinc- 
tion on the "true" (intrinsic) differential magnitude (wt - fih) are 
linearly proportional to X. 

Eventually, we chose the light curve with the smallest 
amount of scatter cr and the lowest level of red noise estimated 
from the jS parameter as defined by Winn et al. (2008). The 
overall S/N of an observation can be quantified by rescaling 
the unbinned photometric cr of the light curve (having a net 
cadence t) on a standard 120-s timescale, that is by calculat- 
ing cri2o = cr Vt/120*. With the only exception of N3, the fi- 
nal IAC80 fight curves have a^o = 0.67-0.95 mmag (Table 1, 
Fig. 2), only sUghtly larger than that achieved by Gibson et al. 
(2008), Tripathi et al. (2010) and Christiansen et al. (201 1), with 
space-based or much larger facilities. N6 shows a large amount 
of red noise of unknown origin, but probably related to color- 
dependent systematics caused by variable atmospheric extinc- 
tion. For the above-mentioned reasons, N3 and N6 were em- 
ployed in the determination of our new ephemeris, but not in 
our TTV analysis (Section 5). 

4.2. UDEM photometry 

Standard dark current subtraction and twiUght sky flat-field di- 
vision process were performed for caUbration on each image of 
the UDEM fight curves (Ul-5). Aperture diff'erential photome- 
try was carried out on the target star and 4-6 comparison stars 
of similar magnitude (\Am\ < 1 .5). The apertures used varied for 
each date due to defocus and weather conditions, but they were 
optimized to minimize the scatter of the resulting fight curves. 
We found that the best results were obtained by averaging the ra- 
tios of WASP-3b to each comparison star. This produced smaller 
scatter than the method of ratioing the target star to the sum of 



all the comparison stars. We estimated the formal error for each 
photometric point as the standard deviation of the ratio to the 
individual comparison stars, divided by the square root of their 
number (error of the mean). 

After normalizing the target star to the comparison stars and 
averaging, some long-term systematics as a function of time 
were found. This is perhaps caused by differential extinction be- 
tween the transit and comparison stars, which generally have 
different and unknown spectral types. This variation was re- 
moved by fitting a Unear airmass-dependent function to the out- 
of-transit basefine of the fight curve. 

5. Data analysis 

5.1. Fitting of fhe transit model 

We chose to analyze all of the new and archival light curves em- 
ploying the same software tools and algorithms. Our goal was 
to get homogeneous estimate of the physical/orbital parameters 
of WASP-3b. In particular we were interested for our TTV anal- 
ysis in estimating the To of each transit and its associated error 
in the most accurate way, avoiding biases due to different tech- 
niques. We avoided estimating Tq through heuristic algorithms 
which assume a symmetric light curve, such that developed by 
Kwee & van Woerden (1956). The main reason is that they do 
not fit for any quantity other than Tq. Even more important, they 
are not robust against outUers, they provide values of Tq known 
to be biased (Kipping 2010) and errors on Tq known to be un- 
derestimated (Pribulla et al. 2012). 

JKTEBOP^ (Southworth et al. 2004) is a code which mod- 
els the light curve of a binary system by assuming both com- 
ponents as biaxial elUpsoids and performing a numerical inte- 
gration in concentric annuli over the surface of each component. 
JKTEBOP version 25 was run to fit a model light curve over our 
data and to derive the four main photometric parameters of the 
transit: the orbital inclination /, the ratio of the fractional radii 
kr = Rp/R*, the sum of the fractional radii = Rp/a + /?*/a 
(Ri, is the stellar radius, Rp the planetary radius, and a the orbital 
semi-major axis), and the mid-transit time T{). We chose to fit 
i, kr, 2r independently for each data set, because in a perturbed 
system i and could change over a long timescale, while kr is 
an important diagnostic of fight curve quafity: when the photo- 
metric aperture is contaminated with flux from neighbors, the 
transit is diluted and k,- becomes smaller. Moreover, an indepen- 
dent fit allows us to highlight correlations between i, kr, 2r and 
to derive more reUable global results, as discussed at the end of 
this Section. 

We set a quadratic law to model the limb darkening (LD) 
efl'ect, naming mi the linear term and U2 the quadratic term: 
Ijj/^o - ^ ~ ''i^l ~ j") ~ ^'2(1 - fi)^ and = cosy, where Iq is 
the surface brightness at the center of the star and y is the angle 
between a line normal to the stellar surface and the line of sight 
of the observer. Southworth (2010), among others, has shown 
that fixing the values of both u i and U2 should be avoided, as it 
could lead to an underestimate of the errors. On the other hand, 
most fight curves have a S/N too low to let both ui and U2 free, 
and the resulting best-fit results can be unphysical. We set the 
quadratic term U2 always fixed to its theoretical value interpo- 
lated from the tables computed by Claret (2000) (BVRJc bands), 
and Claret (2004) (Sloan ugriz), adopting the stellar parameters 
of WASP-3 derived by PoUacco et al. (2008). For all light curves 
from nonstandard photometry, that is unfiltered CCD photom- 
etry (CI -8, LI -3, El) or from wide-band R + V photometry 

^ http : //www . astro . keele . ac . uk/" j kt/codes/ jktebop . html 
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(Gl-2), estimating first-guess LD coefficients is not trivial. As 
for Gl-2, we interpolated ui = 0.24 and U2 = 0.38 from the 
tables by Claret (2000) by taking the average of the values tabu- 
lated for the Johnson-y and Cousin-/? bands, as done by Gibson 
et al. (2008). We did the same for Cl-8, Ll-3, and El, assuming 
that the quantum efficiency of a typical unfiltered CCD usually 
peaks somewhere in between those two bands. 

Then one of the three following procedures was applied: 

1. On the data sets with a high overall S /N and with two or 
more transits gathered with the same instrument and filter 
(Gl-2; Tl-3; T5-T6; Ml-3; M4-5; Cl-8; Ll-3; Nl-5) we 
first fitted a model with free kt, S,, Tq (and ui fixed at 
its theoretical value) to the "best" individual Ught curves to 
get a preliminary estimate of their Tq. For "best" we mean 
complete transits with high S/N: our choice is summarized 
in Table 2, fourth column. Then we phased all those curves 
setting To = 0. The free parameters i, k^, S^, ui were fitted 
again on the stacked light curve, in order to get a high-S/N 
"reference" model of the transit by integrating the informa- 
tion contained in the whole set. We fixed ;, k^, E^, mi to their 
best-fit values, and fitted a model with only Tq as free pa- 
rameter on all the individual transits, including the low-S/N 
or partial ones. 

2. On T4, a high-S/N but single light curve, we carried out one 
simultaneous fit with ;, k^, 2,, Tq, and mi as free parameters. 

3. In all other cases, the data quahty did not allow us to con- 
strain u\ to values with physical meaning, thus u\ was fixed 
to its theoretical value along with U2- Each transit was then 
fitted individually to get kr, E^, and Tq. 

The construction of the reference model from the IAC80 best 
data set (Nl-2, N4-5 light curves) is summarized in (as an ex- 
ample) Fig. 2 It is worth noting that on average the individ- 
ual Ught curves have ctuo = 0.75 mmag, which decreases to 
o"i20 = 0.39 mmag on the stacked data points. The expected 
noise from Poissonian statistics is 0.75/ V4 = 0.375 mmag. By 
taking also into account that our observations do not cover al- 
ways the same orbital phases, we are confident that the level of 
red noise in the IAC80 photometry is very low. 

As the formal errors derived by least squares techniques are 
known to be underestimated in presence of correlated noise, we 
took advantage of two techniques implemented in JKTEBOP to 
estimate realistic errors: a Monte Carlo test (MC) and a boot- 
strapping method based on the cychc permutations of the resid- 
uals (RP or "prayer bead" algorithm, Southworth 2008). The 
errors on all parameters obtained with the RP algorithm are 
on average significantly larger for most of the archival light 
curves, suggesting a non-negligible amount of red noise. We thus 
adopted conservatively the RP results in our subsequent analy- 
sis. Mean values and error bars can be estimated in two different 
ways: 1) as the arithmetic mean of the RP distribution associated 
to its standard error ±cr, and 2) as the median of the RP distribu- 
tion along with its 15.87* (cr_) and 84.13* percentile (o-+). The 
first estimate assumes a Gaussian distribution, while the latter is 
purely empirical: they should match in absence of red noise. 

We adopted as final results the RP median and uncertainties 
cr+, (T_ for every fitted parameter except Tq. The estimated Tq 
will be analyzed in Section 5.2 with periodogram techniques that 
cannot deal with asymmetric error bars, thus for this parameter 
we adopted the RP means with Gaussian errors ±cr. The best-fit 
values of S,, k,, i and ui modeled on each individual data set are 
summarized in Table 2. We show also the theoretical value of 
the linear LD coefficient mi th as interpolated from Claret (2000, 
2004). On all subsets, ui and Mi_th are in agreement within ~ 



Icr. This holds even for Cl-8, Ll-3, and Gl, demonstrating that 
assumptions previously made on nonstandard "clear" or R + V 
photometry are reasonable. 

In an unperturbed system, the transit parameters E^, k^, and 
i are purely geometrical and should not depend on wavelength 
or observing technique. Even if the system is suspected of be- 
ing perturbed, one can check the long-term consistency of E,, ki, 
and i by comparing the best-fit values estimated for independent 
data sets (Table 2). It is then possible to integrate all the ex- 
tracted information to obtain final quantities of higher precision. 
We computed the weighted means of all subset estimate of Er, 
kj, and / Usted in Table 2, obtaining the result shown in the last 
but one line of the same table ({weighted mean)i). As pointed 
out from previous works (e.g., Southworth 2008), these three 
quantities are correlated with each other, as it becomes evident 
by plotting their individual estimate on planes projected from the 
three-dimensional parameter space (E^, k^, i) (Fig. 3). Apart from 
this, the consistency among all measurements is assessed within 
the error bars, and we conclude that no variation of Ei-, ki-, or / 
is detectable over the timescale covered by our data. The data 
point extracted from subset Tl , T3 appears as the only probable 
outlier in two correlation plots out of three. We re-evaluated the 
weighted mean after removing that point ((weighted mean)2 in 
Table 2). The resulting averages are quite similar to those eval- 
uated without removing the outher, but the uncertainty on is 
smaller. We adopt the second mean as final estimate: 



Er = 0.2187 ±0.0098 
kr = 0.1058 ±0.0012 
i = 84°.12±0.82 



(5) 
(6) 
(7) 



More in general, results obtained from subsets Tl-3 (Sloan /) 
and T5-6 (Sloan g) are slightly but significantly in disagreement 
with each other, a fact already noted by Tripathi et al. (2010) and 
then attributed to the presence of residual red noise. 

5.2. TTV analysis 

The best-fit Tq values for each transit, after being uniformly con- 
verted to BJD(TDB), are shown in Table 3 along with their esti- 
mated ± Icr uncertainties (second column). For completeness we 
also tabulated the median value of To estimated from the distri- 
bution of the RP residuals (third column of Table 3) along with 
its 15.87* (cr_) and 84.13* percentile (cr+). 

To check which light curves within our sample are signifi- 
cantly affected by red noise, we consider two diagnostic param- 
eters named jS and E. The first one is defined as in Winn et al. 
(2008): the light curve is averaged over M bins containing A'^ 
unbinned points each, then /3 is calculated as 



<5"N 



M 



(8) 



That is, p is the ratio between the scatter ctn measured on a given 
temporal scale Af = Nt, and the expected noise ctn estimated 
by rescaling the unbinned scatter cr assuming Gaussian statis- 
tics. Ideally, we expect jS ^ 1 for independent and random errors 
(i.e. pure "white noise"), larger values meaning presence of red 
noise at time scales ~ Af. The time scales around Af ^ 25 min 
are the most important ones to our purposes, because they corre- 
spond to the duration of the ingress/egress part of the WASP-3b 
light curve. As shown by Doyle & Deeg (2004), those are the 
parts having the largest information content about To- We chose 
to compute yS on a set of averaging times between Af = 20 and 30 
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IND. TRANSITS RESIDUALS STACKED MODEL 




phase phase phase 



Fig. 2. Construction of tiie best-fit model from the four best light curves of WASP-3b observed at IAC80 (on 201 1 May 7, May 21, Jun 26, Jul 20). 
The ID# of each transit (Nl, N2, N4, N5) matches the corresponding entry in Table 1. Small gray dots represent the data points with the original 
cadence, while blue circles are binned on 120 s intervals. The red line is the best-fit from JKTEBOP (Table 2). Left panel: individual light curves. 
Transits have been offset in magnitude by a multiple of 0.02 for clarity. Middle panel: Residuals from the best-fit model. The reported scatter is 
evaluated on the binned points as the 68.27"' percentile from the median value. Right panel: Stack of all four IAC80 light curves with the best-fit 
model superimposed. The derived parameters are quoted in Table 2 



minutes, then we took their arithmetic mean as a final estimate 
for/3. 

The second diagnostic is E. We mentioned above that the 
distribution of the RP residuals around the best-fit value should 
be symmetric if the noise budget is dominated by white noise. 
Hence a "skewed" distribution could highlight a non-negligible 
amount of red noise. The opposite is not always true: short-term 
systematics (Af < t) do not necessarily lead to skewed RP resid- 
ual distributions. We parametrized this "skewness" with the ratio 
2 between the largest and the smallest error bar cr of a given data 
point: 

E - max{ cr+, cr_)/ min{ cr+, cr_) (9) 

In principle, S ^ 1 for well-behaved transits, and E » 1 for 
transits dominated by long-term systematics. Table 3 lists E for 
all the employed Tq. We found 1.01 < E < 1.63 for the eleven 
TASTE transits (Nl-6, Ul- 5). Instead, a few archival data points 
show unusually large values (e.g., E - 8.07 for 12). We investi- 
gated this issue by comparing the most significant Tq published 
in the literature (Tl-5; Gl-2; Cl-8; Ml-6) with those derived by 
our reanalysis (Fig. 4, third panel from the top). We concluded 
that the vast majority of the published estimate agree within the 
error bars with ours. Notable exceptions are Tl, 12, and G2: they 
all are partial transits, and T2, G2 have also large E. This is ex- 
actly what we expected. It indeed demonstrates that when a light 
curve lacks the off-transit part, its normalization becomes tricky, 
and even a very small difference in the adopted technique can 
lead to significantly different Tq measurements. We emphasize 
this conclusion because many TTV studies employ partial tran- 
sits (some of them even for the most part: among others. Pal 
et al. 2011 and Fulton et al. 2011). While this is fine when es- 
timating the orbital and physical parameters of the planet, we 



demonstrated that partial light curves should be included with 
extreme caution in a TTV analysis. We note that, on average, 
our error bars are larger than the published ones -sometimes by 
a factor of two- confirming our concern that most measurements 
carried out in the past have been published with underestimated 
errors due to neglected red noise. 

We considered two different samples of measurements for 
our TTV analysis. The first ("ALL") includes all the 49 Tq listed 
in Table 3, plus the from the ephemeris Eq. (1) from Pollacco 
et al. (2008). The second sample ("SELECTED") is a high-quality 
subset of 36 values, selected by excluding Eq. (1) (it does not 
correspond to an independent measurement of one single transit) 
and other 13 data points using the following rejection criteria: 

1. partial light curve: data points lacking between the first and 
the last contact; 

2. large scatter: c^^,^) > 1 mmag (0-30 defined as the RMS of the 
residuals after averaging over 30-min bins); 

3. red noise: E > 2 oryS > 1.2; 

4. presence of systematics during or close to the transit 
ingress/egress, as determined by visual inspection of the 
residuals. 

The transits excluded from SELECTED are listed in Table 3 along 
with the reason for the exclusion. 

The first step to plot an (9 - C diagram is to calculate a "ref- 
erence" linear ephemeris to predict To at any given epoch. We 
set the new zero epoch at N2, i.e. our most accurate light curve. 
The ALL sample was employed to fit a linear model by ordinary 
weighted least squares, obtaining 

ro(BJDTDB) = 2 455 702.57993 ±0.00017 + 

+N ■ 1.8468349 ± 0.0000004 (10) 
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Fig. 3. Geometrical parameters of WASP-3b estimated from individual data (sub-)sets (Table 2) plotted as black labeled circles in their two- 
dimensional parameter space, to highlight the sizeable correlation between Ei-, fe,., and /. The blue dashed line in each plot is a weighted linear fit 
of all points; the blue triangles and the associated error ellipse correspond to the weighted mean ± Icr computed on the marginal distribution of E,., 
k^, and i (last but one row of Table 2). The red continuous line and squares are computed in the same way as the blue ones, but after removing the 
outlier labeled "3" from the set (Tl , T3; last row of Table 2). 

Table 2. Orbital/physical parameters of WASP-3b estimated from individual data (sub-)sets. 



reference paper 


Ni 


telescope 


ID 




band 


I,. 


kr 






III 


f l.th 


Gibson et al. (2008) 


1 


LT-2.0m 


Gl 




(R+V) 


"■■^■^^^-0.0064 





1 Ofi«+0.0008 


83° 44+<'<i6 




0.28 


Tripathi et al. (2010) 


2 


FLWO-1.2m 


T5- 


-6 


Sloan g 


2266+° '"'™ 





IOSA+0.0014 


83". 25!°;?^ 


0.380:SJ;»55 


0.42 


Tripathi et al. (2010) 


2 


FLWO-1.2m 


T1,T3 


Sloan / 


2099+° '"'^^ 





1 n 1 c+0.0009 
^"^-^-0.0009 


85°.04:«;^^ 


n 249+° "'*° 


0.18 


Tripathi et al. (2010) 


1 


UH-2.2m 


T4 




Sloan z 


2164+° "™ 





1065+°°°°^ 


84°.56!«;^^ 


141+°"5<^ 

"■^^^-0.054 


0.14 


Maciejewski et al. (2010) 


3 


Rohzen-0.6m 


Ml- 


-3 


R 


2162+°°"^ 





lQgl+0.0014 
^""^-0.0021 


84° .42:1 -J J 


221+0130 


0.24 


Maciejewski et al. (2010) 


3 


Jena-0.9m 


M4- 


-6 


R 


2195+° °'^^ 





1 nc2+0.0014 
^"-^ -0.0013 


83°.89:;i^ 


095+° i^'* 

"■"='-'-0.129 


0.24 


Christiansen et al. (201 1) 


6 


HRI@EP0X1 


Cl- 


-4,C7-8 


clear 


A 22 10+0.0049 





1061+0.0008 
^ -0.0008 


83°.87:«;=i 


270+° "*'' 

-0.054 


0.28 


Littlefield (2011) 


3 


SC-U" 


Ll- 


-3 


clear 


Q 2361+0.0205 





1 nS 1 +0.0033 
-0.0040 


82°.90!i|-] 


455+0.158 

>J-^-'-'-0.183 


0.28 


this work 


4 


1AC80 


Nl- 


-2,N4-5 


R 


n ')n9'!+o oo40 





1046+0.0006 

^"^"-0.0006 


85".24!«;=|; 


247+° "^'' 

"■^^'-0.028 


0.24 


(weighted mean) i 












n 21 7'?+0 0095 
O-^i '-'-0.0095 





1 nS'?+0''019 


84°.24tg 






(weighted mean)2 












n Tl 97+0.009S 





1 o^e+0.0012 


84M2t- 







Notes. The columns give: the reference paper, the number of light curves fitted, the telescope, the ID code of the transits, the filter employed, the 
fitted sum of the fractional radii (E, = R^la + Rf,l a), the fitted ratio of the fractional radii (fc, = R^jR^), the fitted inclination i, the fitted linear limb 
darkening (LD) coefficient , and the theoretical linear LD coefficient wi ,i, interpolated from the Claret (2000, 2004) tables. The quoted error bars 
are derived from the JKTEBOP RP algorithm. The last two rows show the weighted means of all the previous individual estimate (1) and of all of 
the previous with the exception of the data set Tl , T3 (see text for details). 
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The uncertainties have been evaluated from the covariance 
matrix of the fit, and were both rescaled by to take into 
account the real dispersion of the data points around our best-fit 
ephemeris. In Fig. 4 are plotted the 0-C diagrams for ALL (first 
panel from top) and SELECTED samples (second panel, with a 
smaller baseline). In both diagrams, the reduced = 3.7-4.5 
suggests that the measurements are not in fuU agreement with 
the linear ephemeris in Eq. (10). Yet, there is no evident periodic 
pattern in our diagrams. We investigated the possibility that this 
statistically significant scatter is caused by a genuine TTV, either 
the TTV claimed by Maciejewski et al. (2010) or a different one. 

We exploited two algorithms to search for periodic signals: 
the Generalised Lomb-Scargle periodogram (GLS; Zechmeister 
& Kiirster 2009) and the "Fast x^'' algorithm (F;^f^; Palmer 
2009). Both these techniques are able to deal with irregularly 
sampled data with nonuniform weights, and minimize aUasing 
eflfects due to the window function. In addition to this, Fx^ can 
also search for an arbitrary number of harmonics. 

We searched for periodic signals with GLS in the period 
range P - 4-1000 d, the lower limit being imposed by the 
Nyquist sampling criterion to avoid aliasing (Home & Baliunas 
1986). The resulting periodograms for both samples are plotted 
in the fourth panel of Fig. 4. In neither case a prominent peak is 
visible. The ALL and SELECTED periodograms are quite similar, 
and their highest peaks stand at P{k) ^ 6.41 d and P(S) ^ 1 1.29 
d, respectively. The formal false alarm probability (FAP) as de- 
fined by Zechmeister & Kurster (2009) is 0.023 (2.28cr) for the 
ALL power peak and 0.076 (1.77 cr) for the SELECTED peak, 
i. e. only marginally significant. However, the formal FAP is de- 
rived under the assumption of pure Gaussian noise, which is not 
our case. To take into account the intrinsic dispersion of our data, 
we investigated whether these peaks are statistically significant 
or not with a resamphng algorithm. We generated 10000 syn- 
thetic 0-C diagrams with the same temporal coordinates of 
actual data points, by randomly scrambling the 0-C values 
at each generation. A GLS periodogram was then evaluated on 
each of them with the same settings of that applied on real data. 
The power of the highest peak found in the real data for ALL and 
SELECTED samples Ues respectively at the 12* (-1.17o-) and 
38* percentile (-0.3 Icr) of the distribution of the maximum- 
power peaks in the synthetic, randomly-permutated diagrams. 
We conclude that neither peaks can be considered as statisti- 
cally significant. In particular, the P(M) ^ 127 days periodic- 
ity claimed by Maciejewski et al. (2010) is not consistent with 
our data. Instead of a peak, the periodogram range where the 
f (M) peak is expected is characterized by an extremely low GLS 
power (fourth panel of Fig. 4, blue region). On the other hand, 
tests on synthetic 0-C diagrams having the same sampling 
and noise properties of our sample demonstrate that aP = f (M), 
A(0-C) = 0.0014 days signal corresponding to the Maciejewski 
et al. (2010) claim would be easily detectable from our data. 

We set Fx^ to search for periodicities with one or two har- 
monics in the same frequency range. Results from both samples 
are quite similar to those obtained above with GLS, with non- 
significant power peaks at periods very close to those previously 
found in the GLS periodogram. 

We carried out our analysis also on a subset of data points 
corresponding to those analyzed by Maciejewski et al. (2010), 
by employing the same tools used for our full data set (Fig. 5). 
When plotted adopting the same frequency limits, the result- 
ing periodogram (and the peak corresponding to the maximum 
power) perfectly matches those published by Maciejewski et al. 
(2010) (upper right panel of Fig. 5). On the other hand, the up- 
per limit in frequency set by Maciejewski et al. (2010) (0.038 



P ') is too low according to the Nyquist criterion. In the range 
0.038 < y < 0.5 P~^ many other maxima are as high as the 
0.0145 P~^ peak (lower panel of Fig. 5). Following a statistical 
test similar to the one above described, we discard that peak as 
not significant, being due to small-sample statistics. 



6. Discussion and conclusions 

In the present study, we analyzed eleven unpublished light 
curves of WASP-3b and re-analyzed other thirty-eight archival 
hght curves, all of them with the same software tools and proce- 
dures. We derived improved orbital and physical parameters for 
this target (Table 2), and computed a refined ephemeris (Eq. 10). 
All individual measurements of the central instant To have been 
compared with the new ephemeris to search for changes in the 
orbital period P of the transiting planet. We concluded that avail- 
able observations of WASP-3b, spanning more than four years, 
are not consistent with a linear ephemeris (x^ ^ 4). A possible 
explanation for this scatter is the presence of a perturbing body 
in the WASP-3 planetary system. 

It is known that the impact of red noise on high-precision 
transit photometry is still not fully understood. F*revious claims 
of TTVs have been disproved on this basis (Southworth et al. 
2012; Fulton et al. 2011). Could the observed scatter in the 
0-C diagram of WASP-3b be explained in terms of underes- 
timated observational errors or caUbration issues? The absolute 
time calibration of each archival light curve cannot be indepen- 
dently checked. In principle one should trust the authors about 
that. However, we point out two main clues supporting the auto- 
consistence of the overall data, and hence the hypothesis of a 
genuine TTV: 

- on three dififerent epochs (A^ - 444, 486, 653 following the 
PoUacco et al. 2008 ephemeris) multiple observations of the 
same transit are available. As they were carried out by differ- 
ent authors at difTerent facilities, they should represent inde- 
pendent measurements of the same quantity. All these data 
points (marked with a star in the first coluimi of Table 3) 
agree with each other within their Icr error bars, suggesting 
that the uncertainties on Tq are correctly evaluated by our 
pipeline; 

- some anomalous patterns in the O - C diagrams are con- 
firmed by several different data sets. For instance, nearly 
all points gathered in 2009 within the range = 440-510 
lie ahead of the Tq predicted by our baseline ephemeris 
{O - C < 0). The only exception is U4, which essentially 
lies at O - C ~ within its error bar. The weighted mean 
of these twelve measurements from eight different authors is 
O-C = -0.00118 ±0.00016 days. This implies a 7.2cr devi- 
ation from a constant orbital period. These patterns can also 
be detected in high-precision data subsets, such as our ones. 
Among the four best IAC80 transits, the first two (Nl-2) are 
delayed by 4.3cr (i.e. 0-C = 48 ±11 s) compared with 
the prediction, while the second two (N4- 5) are ahead of the 
ephemeris by 2.8cr (O - C = 35 ± 13 s). 

The available data thus strongly suggest an intrinsic devi- 
ation of the actual transit times from the linear ephemeris ex- 
pected from a Keplerian two-body transiting system. On the 
other hand, our analysis rules out the periodic TTV claimed by 
Maciejewski et al. (2010), and failed even at detecting any sig- 
nificant periodicity in the updated O-C diagram. Most TTV 
studies carried out in the past searched for a periodic signal, but 
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Fig. 4. First panel from the top: O - C diagram for all data points tabulated in Table 3. Second panel: same as above, for points selected as in the 
last column of Table 3. Third panel: Comparison between the original Tq published by the respective authors (white circles with gray error bars; 
Table 1) and as re-estimated in this work (red triangles and bars), for the subset of high-precision light curves identified in the horizontal axis. 
Fourth panel GLS periodogram for the complete sample (red line, highest peak at A) and the selected sample (black line, highest peak at 5). The 
periodicity claimed by Maciejewski et al. (2010) is marked with the M label. 
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Fig. 5. Period analysis carried out on the same set of O - C data points analyzed by Maciejewski et al. (2010). Upper left panel: O — C diagram 
for the selected points. Upper right panel: GLS periodogram as a function of frequency y, adopting the same plotting limits of Maciejewski et al. 
(2010). Bottom panel: GLS periodogram as a function of period, adopting wider limits on frequency according to the Nyquist criterion (see text 
for details). The red vertical line marks the peak claimed by Maciejewski et al. (2010). 



here we are dealing with a more complex and non-periodic phe- 
nomenon, at least at the time scales we sampled. What kind of 
dynamical system can induce such a perturbation? 

Veras et al. (2011) demonstrated that some orbital config- 
urations, especially close to (but not exactly in) mean-motion 
resonances, can induce quasi-periodic or even chaotic TTVs. In 
other non-exotic configurations, the periodicity would manifest 
itself only at time scales > 10 yr (Veras et al. 201 1). Also when 
more than one perturber are present and their orbital periods are 
not commensurable, as in the case of our inner Solar System, the 
resulting TTV would be in general aperiodic (Holman & Murray 
2005). 

If WASP-3b belongs to one of the above-mentioned cases, 
careful dynamical modeling and additional follow-up is required 
to confirm the hypothesis and to constrain the mass and period of 
the possible perturber(s). Photometric monitoring is still ongo- 
ing within the TASTE project, and high-precision RV measure- 
ments are foreseen with HARPS-N. As stressed out by Meschiari 
& Laughlin (2010) and Payne & Ford (201 1), photometric TTVs 
and RVs are highly complementary in breaking the degeneracies 
that are common in the inverse dynamical problem. 
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Appendix A: A new eclipsing variable 

A star in the WASP-3 field (UCAC3 a = 18*" 34™ 07\36, 5 = 
+35°38'59".6, epoch 2000.0; Zacharias et al. 2010), initially 
chosen as reference star, was found to be variable and excluded 
from the reference list. A complete light curve of this variable 
was assembled by registering the magnitudes on the Nl-2 and 
N4-N6 frames on a common zero point (Fig. A.l, upper panel). 
When folded on the P ^ 0.3524 days peak detected in the Lomb 
& Scargle periodogram, the binned curve shows a periodical 
pattern with equal maxima and slightly different minima, typi- 
cal of contact eclipsing binaries (W UMa-type; Fig. A.l, lower 
panel). We derived the following ephemeris, setting <l> = at the 
phase of primary minimum and estimating uncertainties through 
a bootstrapping algorithm: 

JoCBJOtdb) = 2455 776.48932 ±0.00024 + 

+ 0-0.3523683 ±0.0000018 (A.l) 

This variable star appears to be unpublished, and we sub- 
mitted it to the International Variable Star Index (identifier: 
VSXJ183407.3+353859). Colors from catalog magnitudes: V = 
15.63, R = 14.99 (NOMAD; Zacharias et al. 2004), J = 14.64, 
H = 14.34, = 14.21 (2MASS; Skrutskie et al. 2006), 
and proper motions /i„cos<5 = -32 mas/yr, ^i^ - 19 mas/yr 
(UCAC3) suggest that this object could be a binary with both 
components of late-G spectral type. 



Nl \:i N4 N5 N6 




phase (P = 0, 352368 d) 



Fig. A.l. Light curve of a previously unreported R ~ \5 variable star 
in the WASP-3 field, classified as a contact eclipsing binary (see text 
for details). Top panel: unbinned data points folded on the best-fit pe- 
riod P = 0.353626 days. Different nights are coded in different colors. 
Bottom panel: same as above, binned on 0.02 intervals in phase. 
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Table 1. Summary of the light curves of WASP-3b analyzed in this work. 



N 


ID 


evening 


telescope 


band 






c 




reference paper 


notes 






date 








fs) 


( mTnaff^ 


V^ixuixag/ 






250 


Gl 


2008/05/18 


LT-2.0m 




4799 


3.0 


3.97 


0.63 


Gibson et al. (2008) 




309 


G2 


2008/09/04 


LT-2.0m 


(R + V) 


3900 


3.0 


3.81 


0.61 


Gibson et al. (2008) 


partial 


248 


Tl 


2008/05/15 


FLWO-1.2m 


Sloan / 


403 


33.8 


1.38 


0.73 


Tripathi et al. (2010) 


partial 


262 


T2 


2008/06/10 


FLWO- 1 .2m 


Sloan i 


269 


33.6 


1.16 


0.61 


Tripathi et al. (2010) 


partial 


268 


T3 


2008/06/21 


FLWO-1.2m 


Sloan i 


522 


34.2 


1.59 


0.85 


Tripathi et al. (2010) 




280 


T4 


2008/07/13 


UH-2.2m 


Sloan z 


276 


67.5 


0.81 


0.61 


Tripathi et al. (2010) 




444 


T5 


2009/05/12 


FLWO- 1 .2m 


Sloan ^ 


224 


78.3 


1.17 


0.95 


Tripathi et al. (2010) 




451 


T6 


2009/05/25 


FLWO- 1 .2m 


Sloan ^ 


238 


74.5 


1.08 


0.85 


Tripathi et al. (2010) 




486 


Ml 


2009/07/28 


Rohzen-O 6m 




283 


49.3 


2.06 


1.32 


Maciejewski et al. (2010) 




499 


H2 


2009/08/2 1 


Rohzen-O 6m 




200 


72.2 


2.36 


1.83 


Maciejewski et al. (2010) 




506 


H3 


2009/09/03 


Rohzen-O 6m 




297 


45.8 


2.87 


1.78 


Maciejewski et al. (2010) 




519 


M4 


2009/09/27 


Jena-0 9m 




259 


59.6 


2.72 


1.92 


Maciejewski et al. (2010) 




539 


M5 


2009/1 1/03 


Jena-0 9m 




258 


51.6 


2.41 


1.58 


Maciejewski et al. (2010) 




629 


M6 


2010/04/18 


Jena-0 9m 




206 


78.0 


1.50 


1.21 


Maciejewski et al. (2010) 




486 


Dl 


2009/07/28 


OAVt1A-0 25m 


f( 


353 


46.7 


2.30 


1.43 


I../ullldKi>U CI til. \Z.\JL\J1 




290 


CI 


2008/07/3 1 


HRIOFPOXI 




1 157 


53.2 


0.89 


0.59 


Christiansen pt al ^901 H 




291 


C2 


2008/08/02 


HRT®FPOXT 

X XXVI. x_rfXr v-«vx 


clear 


958 


51.7 


1.01 


0.66 


Christianspn pt al f'901 H 




292 


C3 


2008/08/04 


HRI@EPOXI 


clear 


985 


51.3 


1.03 


0.67 


Christiansen et al ('2011't 




293 


C4 


2008/08/06 


HRI@EPOXI 


clear 


1138 


51.3 


0.98 


0.64 


Christiansen et al ("201 




294 


C5 


2008/08/08 


HRI@EPOXI 


clear 


281 


52.2 


0.97 


0.64 


Christiansen et al ("201 Ti 


nartial 

L'tU UtU 


296 


C6 


2008/08/1 2 


HRI@EPOXI 


clear 


954 


54.6 


0.93 


0.63 


Christiansen et al nO1 1^ 


nartial 

L/tU UtU 


297 


C7 


2008/08/13 


HRI@EPOXI 


clear 


640 


51.5 


1.02 


0.67 


Christiansen et al. (2011) 




298 


C8 


2008/08/15 


HRI@EPOXI 


clear 


852 


51.5 


1.07 


0.70 


Christiansen et al. (2011) 




654 


LI 


201 0/06/04 


SC-11" 




218 


76.8 


2.68 


2.14 


Littlefield (2011) 




686 


L2 


2010/08/02 


SC-ll" 


clear 


192 


113. 


2.38 


2.31 


Littlefield (2011) 




706 


L3 


2010/09/08 


SC-11" 


clear 


192 


65.5 


2.70 


1.99 


Littlefield (2011) 




488 


Zl 


2009/08/01 


Weihai- Im 


V 


323 


43.4 


7.35 


4.42 


Zhang et al. (2011) 




444 


SI 


2009/05/12 


VCT-0.5m 


Sloan 7, 


177 


94.7 


2.25 


1.99 


Sada et al. (2012) 




653 


S2 


2010/06/02 


KPNO-2.0m 


J 


250 


71.8 


1.50 


1.16 


Sada etal. (2012) 




653 


S3 


2010/06/02 


VCT-0.5m 


Sloan z 


186 


91.1 


2.23 


1.94 


Sada et al. (2012) 




842 


S4 


2011/05/17 


VCT-0.5m 


Sloan z 


217 


87.4 


3.52 


3.00 


Sada et al. (2012) 




194 


El 


2009/04/22 


Newton-0.2m 


clear 


194 


74.2 


3.36 


2.64 


ETD obs. Trnka 




446 


E2 


2009/05/16 


SC-12" 


/ 


181 


104. 


3.18 


2.96 


ETD obs. Gregorio 




653 


E3 


2010/06/02 


SC-12" 


R 


177 


72.7 


2.90 


2.26 


ETD obs. Shadik & Patrick 




666 


E4 


2010/06/26 


Newton-0.3m 


B 


266 


61.6 


3.69 


2.64 


ETD obs. Garlitz 




699 


E5 


2010/08/26 


RC-12.5" 


R 


119 


124. 


2.62 


2.67 


ETD obs. Hose 




838 


E6 


2011/05/09 


Monteboo-0.6m 


R 


393 


35.3 


2.55 


1.38 


ETD obs. Janov 




849 


E7 


2011/05/30 


RC-12" 


V 


297 


51.2 


5.65 


3.69 


ETD obs Dvorak 




837 


Ul 


2008/09/03 


UDEM-0.36m 


J 


467 


33.7 


4.27 


2.26 


this work 




844 


U2 


2009/06/05 


UDEM-0.36m 


I 


566 


34.1 


3.78 


2.01 


this work 




864 


U3 


2009/07/25 


UDEM-0.36m 


Sloan z 


383 


45.1 


4.38 


2.68 


this work 




877 


U4 


2009/08/05 


UDEM-0.36m 


Sloan z 


371 


44.2 


5.03 


3.05 


this work 


lowS/N 


884 


U5 


2010/08/15 


UDEM-0.36m 


7 


471 


35.0 


3.31 


1.79 


this work 




837 


Nl 


2011/05/07 


IAC-0.8m 


R 


678 


24.0 


1.78 


0.79 


this work 




844 


N2 


2011/05/21 


IAC-0.8m 


R 


742 


23.1 


1.53 


0.67 


this work 




851 


N3 


2011/06/02 


IAC-0.8m 


R 


412 


38.1 


2.93 


1.65 


this work 


lowS/N 


864 


N4 


2011/06/26 


IAC-0.8m 


R 


566 


25.1 


1.53 


0.70 


this work 




877 


N5 


2011/07/20 


LAC-0.8m 


R 


785 


20.1 


2.12 


0.87 


this work 




884 


N6 


2011/08/02 


IAC-0.8111 


R 


694 


19.9 


2.34 


0.95 


Ihis work 


red noise 



Notes. The columns give: the transit epoch N assuming T = Tq+NP and the original ephemeris (Eq. 1) from PoUacco et al. (2008), the ID code, the 
"evening date" of the observation, the filter employed, the number of unbinned data points, the average net cadence t in seconds, the photom etric 
scatter cr measured as the 68.27* percentile of the residuals from the best-fit model, the normalized photometric scatter cri2o = crVT/120% the 
reference paper or database, and comments. 
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Table 3. Central instants of WASP-3b transits estimated from all the individual light curves. 



N 


io (mean), J3JU(1UJ5) 


To (med), BJD(TDB) 




O — C (s) 




ID 


telescope 


selected 


(0) 


2454143. 85107 ± 0.00040 




1.00 


A AAA 1 A 

-0.(XK)19 


-0.50 






no (baseline ephemeris) 


248 


2454601. 800/1 ±0.00026 


'>AZA£.r\\ o<;<;-7i+l) 00025 

245460 1 .866 / -J_n qqq^^ 


1. 12 


A AAA1 O 

0.(X)038 


1.48 


Tl 


T^T A17/^ 1 

rLWU-1.2m 


no (partial) 


250 


2454605.56042 ± O.OOOJO 


2454605 56049+"-''''''^'* 


2.00 


A AAA A n 

0.(XXJ42 


1.42 


Gl 


LT-2.0m 


yes 


262 


245462/. /2j1 / ± 0.00068 


2454627.72324_„ 5 


8.07 


0.001 15 


1.70 


T2 


rLW U- 1 .2m 


no (partial) 


268 


24546J5. 80399 ± 0.000J4 


A Z A /L'~> O OA.1AA-H0 OOO^l) 

245463 8 . 80400_g ^^p^y 


1.27 


A AAAA£ 

o.(xxjy6 


2.85 


T3 


T7T A17/^ 1 

rLW U- 1 .2m 


yes 


280 


^AZA/i/Lf\r\/LAnc\ f\ f\r\f\ 1 c 

2454660.964/9 ± 0.00015 


A z A £ £r\ r\£ A or\+\i 00(3 15 

2454660. 96480_yyyyJ^ 


1.00 


A AAAO C 

-0.(XXJ25 


-1.68 


T4 


UH-2.2m 


yes 


290 


"1 /I c /I £.nc\ /I 1 1 ri f\f\r\ a o 
24546 /y.4ii 18 ± 0.00042 


2454679 433 1 1 +o-ooo55 

'^ACAiio^ mn/CQ+O 00038 

2454681.27963_p„„„i 


1.90 


A AAA1 1 

-0.00021 


-0.50 


CI 


rU<I<S't,rUAl 


yes 


291 


245468 1.27967 ± 0.00034 


1.23 


A AAACC 

-0.00055 


—1.63 


C2 


T TT> T I~" TI/'A'V 1 

HK1@EPUa1 


yes 


292 


2454oo3. 12798 ± 0.00049 


2454683.12800!°;™; 

O/lC/liiQ/l ATCT 1 +0.00042 

2454oo4.9/52i_Q 


1.09 


A AAAA1 

0.00092 


1.88 


C3 


T Tr> T iOv T7"r*/^VT 


yes 


293 


2454oo4.9/j24 ± U.UUU4U 


1.11 


A AA1 1A 

0.00134 


3.36 


C4 


tTD T /SlCTli^VT 

liKI (S'c.r'UAl 


yes 


294 


24j4ooo.o2144 ± U.UUU/4 


O/lC/liCQ^ Q01 /I A+O.OOiOz 

245468b.82140_p „p|)4. 


2.17 


A A A ATI 
0.000 / 1 


0.96 


C5 


tTD T /SlCTli^VT 

rLKI (SPcrUAl 


no (partial) 


296 


2454o90.5i4/5 ± 0.UUU59 


O/lC/liiAA CI /lOA+0. 00048 

2454690.51480 

l/lC/liiAl '3*;iTi+0. 00047 

2454692.36173_„ 


1.46 


A AAA^C 

0.00035 


0.59 


C6 


UTJ T /SlCTli^VT 

rlKI (SPcFUAl 


no (partial) 


297 


2454o92.3olDo ± 0.00056 


1.36 


A AAA/1 A 

0.00044 


0.79 


C7 


HKI (cp trUAl 


yes 


298 


2454694.20 //6 ± U.UUUoJ 


i/iC/i(CO/i lA'/^n+0. 00090 
2434694.20 / /O ^ 01)1)7:^ 


1.20 


A AAA1 1 

— U.(XX13 1 


—0.37 


C8 


HKlvff iil:'(jAl 


yes 


308 


2434/12.6/641 ± U.(XXJ64 


2454712 6 763 7 


1.06 


A AAAAA 

— O.(XXXa) 


—0.01 


Ul 


UUt:,JVl-0.3om 


yes 


309 


'lA^AI \ A ^nif^Q _i_ n f\f\f\A 1 

Z4j4/ 14.JZJDfi ± U.UUU41 


2454714.52358!; « 

"i/iC/in/iQ cio/i n+0. 00045 

2454943.53249_„„„„,4 


3.15 


A AAA/1 1 
U.00U4Z 


1 .04 


G2 


LT-2.0m 


no (partial) 


433 


'^AZAC\A1 zi"iAr\ 1 f\ r\r\f\^r\ 
2454943. 3J24U ± U.(XXJ6U 


1.64 


0.00161 


2.70 


El 


Newton-0.2m 


yes 


444 * 


2434963.84430 + U.(XX)81 


2454963. 84453!°°°°g 
2454963 84541+°-°°"™ 


1.16 


A AA 1 A C 

— 0.(X)146 


— 1.81 


T5 


C7T \li7/^ 1 1.™ 

rLWU-1.2m 


yes 


444 * 


2434963.8432/ +0.1X)1I8 


1.31 


A AAA/; A 

— u.o(XJ6y 


—0.59 


SI 


1 Tf^'y A C™ 

VL.l-0.3m 


yes 


446 


2454967.53905 ± 0.00070 


^ A Z AC\/1.^ COA1^-HO OOO6I 

2454967 .5391 2!n;QQ™ 


1.30 


A AAAC O 

-0.00058 


-0.84 


E2 


SC-12" 


yes 


451 


2454976.77284 ± 0.00030 


1/ICy)A"7Z: T71OO+0 000l2 

2454976.77283_„„„„28 


1.14 


—0.00097 


—3.23 


T6 


TTT \"I It^ 1 1«« 

rLWU-1.2m 


yes 


457 


2454987.85256 ± 0.00093 


2454987.85267+°°°°?' 


1.63 


-0.00225 


-2.43 


U2 


UDEM-0.36m 


yes 


484 


Z4550J/./io/o ± O.OOOiSo 


2455037.71877!; 


1.01 


A AAACO 

— 0.00058 


—0.68 


U3 


U JJcM-0. 3om 


yes 


486 * 


2455041.41172 ± 0.00035 


1/1 CCA/1 1 d 1 1 Tl+O 00036 

2455041.41172 „ „„p3, 

l/t CCA/1 1 /I 11/1 *C+0. 00053 

2455041.41246_ 

1 /I CCA/1 C 1 AC C "7 +0.00 1 12 

2455045.10557 „„„,„„ 


1.03 


A A A 1 1 1 

—0.00131 


—3.75 


Dl 


(JAVaA-0.25m 


yes 


4oO * 


Z433U4i.41Z33 ± U.UUU3o 


1 .19 


A AAA/1 Q 
— U.00U4O 


—0.83 


Ml 


Rohzen-0.6m 


yes 


488 


tAZZf^AZ 1 r\z^z 1 A AA 1 rii; 

2433043.10363 ± 0.(X)1U6 


1.12 


A A A 1 AC 

— U.(X)1U3 


—0.99 


Zl 


Weihai-lm 


no (low o/JN) 


490 


1/1 CCA/1 Q QAA^iC _i_ n AA 1 ClI 
Z43304fS.oOUD3 ± U.OOlU/ 


1/1CCA/1Q QAA^^A+0. 00104 
Z43304fS.oOUD6 qiiqiiq 

2455065.42026!™°;™ 
2455078 34799+° °°'-'-' 


1 .06 


A AAA1"7 
U.OOUZ / 


0.26 


U4 


UUilJVl-0.36m 


no (red noise) 


499 


1/1 CCA/^C /I imi _i_ A AAAIA 
Z433063.4ZUZ3 ± U.00U36 


1.12 


—0.00165 


—4.60 


112 


Rohzen-0.6m 


yes 


506 


1/1 CCATQ 1/1 QflO _i_ A AA 1 1 A 
Z4330/o.34fSUy ± U.00il4 


1.39 


—0.00163 


— 1.44 


113 


Rohzen-0.6m 


yes 


519 


1/1 C CI AO QAAin _i_ A AAAO/I 
Z43310Z.36U3U ± U.00Uo4 


2455102 16037+°°°°''' 


1.75 


0.00171 


2.04 


M4 


Jena-0.9m 


yes 


539 
ozy 


T/1CC1 IQ lOTCI J. A AAA"71 

Z4j313y.Zy /j3 ± U.OOU/3 
i/iccanc ci 1 1 T J. A aaac/^ 


2455139.29719!°°°i;i 

1/1 cc^nc cm 0+0:00049 
2455305.51119 nooo.^ 
i/iccri/in o-aoi (;+o.00085 
Z43j34y. 63316 nnQ096 


3.26 
1.29 


A AA11/1 
U.00ZZ4 

A Art AT /I 
U.UUU/4 


3.08 

1 lA 
1.34 


M5 
H6 


Jena-0.9m 
Jena-0.9m 


yes 
yes 


Odd * 


'^A«1AQ Q1'if\A _i_ A AAAQA 

z43jo4y.ojjUo ± u.ouuyu 


1 1 1 

1.13 


A AAl AC\ 
— U.UU14U 


1 c/; 


S3 


\//~'T' A Cm 


yes 


653 * 


'^AZCIAn OllCif^ _i_ A AAA<Cn 

2455349.00390 ± 0.00069 


2455349.83384:1™ 


1.00 


A AAAC/: 
— U.UUU50 


—0.81 


S2 


JsJ:'JNU-2.Um 


yes 


653 * 


'^ACO Af\ O'y Al A ■ A AAAC A 

2455349.83434 ± 0.00054 


i^ACcy Af\ 01 /fjA+O 00061 

2455349.83430!n™j;7 


1.30 


A AAA1 1 

—0.00012 


—0.22 


E3 


SC-12" 


yes 


654 


1/tCC1Cl A t r\ r A A A ICC 

2455351.68410 ± 0.00155 


1/1CC1C1 ^c^An+O 00 [ 5 1 

2455351.68399!";""',' 


1.05 


A AAO OA 

0.00280 


1.81 


LI 


SC-ll" 


no (red noise) 


666 


l/ICC^^O O A'^OCl A AAA A Z 

24553 /3. 84289 ± 0.(XX)45 


A z znn't Ovi'^oc +0.00048 

2455373.84285 , 

^ A Z Z A lA 'TOAOT+O 00148 

2455410.78083!";""™ 

'i.icC/in TAOAi+O 00050 

2455423.70891!";"""" 

"\ A Z Z An A ^OT7C+0 00041 

2455434.78775 

Z43344/./10Z0_„yy^ 


1.12 


A AAA A ^ 

-0.(XX)42 


-0.94 


E4 


Newton-0.3m 


yes 


686 


I/ICC/IIATOATT A A A 1 T T 

2455410. /80/3 ± 0.(X)1 /3 


1.34 


A AAAT 1 

0.000/1 


0.42 


L2 


SC-1 1" 


no (red noise) 


693 


2455423. /0889 ± 0.(XXj48 


1.06 


A AA 1 AO 

0.00103 


2.15 


U5 


UUbM-0.36m 


yes 


699 


^ A Z Z A'~t A 10^^^ A AAAC 1 

2455434. 1>SI 11 ± 0.(XX)51 


1.46 


A A A 11/1 

-0.(X)1 14 


-2.25 


E5 


RC-12.5" 


yes 


706 




3.61 


-D 0007 1 


-0.61 


L3 


SC-1 1" 


HA ( rpri noic^'^ 


837 


2455689.65263 ± 0.00015 


2455689.65263+°™ 


1.14 


0.00054 


3.66 


Nl 


lAC-0.8m 


yes 


838 


2455691.49899 ± 0.00064 


245569 1.49899+§™ J 
2455698.88463+g™ 


1.27 


0.00007 


0.11 


E6 


Monteboo-0.6m 


yes 


842 


2455698.88476 ± 0.00160 


1.18 


-0.00149 


-0.94 


S4 


VCT-0.5m 


yes 


844 


2455702.58052 ± 0.00028 


2455702.58051+ J™ 


1.04 


0.00059 


2.12 


N2 


IAC-0.8m 


yes 


849 


2455711.81615 ±0.00107 


2455711.81617« 


1.28 


0.00204 


1.92 


E7 


RC-12" 


no Oow S/N) 


851 


2455715.50882 ± 0.00057 


2455715.50853+j!:SS 


1.24 


0.00105 


1.84 


N3 


IAC-0.8m 


no Oow S/N) 


864 


2455739.51620 ±0.00017 


2455739.51620+;ii;;;|^ 


1.20 


-0.00042 


-2.49 


N4 


IAC-0.8m 


yes 


877 


2455763.52511 + 0.00031 


2455763.52511 ' 


1.07 


-0.00036 


-1.19 


N5 


lAC-0.8m 


yes 


884 


245577().4.S|y2 + O.OOOS') 


2455776.45192 


1.0') 


-O.OOUO 


-1.58 


N6 


lAC-O.Siii 


no U'clI iKiisc) 



Notes. The columns give: the transit epoch N assuming J = To + NP and the original ephemeris from PoUacco et al. (2008), the best-fit value for 
the central instant Jq of the transit and the associated cr, the median value of the distribution of Tq from the Residual Permutation (RP) algorithm 
and the associated errors cr+ and it_, the "skew" parameter S - max( cr+, it_)/ min((T+, cr_), the O - C according to the new ephemeris (Eq. 10) in 
seconds, the O - C in units of cr, the ID code of the light curve, the telescope employed, and comments about whether and why the data point is 
excluded from the "selected" sample. Simultaneous transits are marked with a star in the first column. 



